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Abstract 

A systematic study of the temperature and magnetic field dependence of the longitudinal and 
transverse resistivities of epitaxial thin films of magnetite (Fe304) is reported. The anisotropic 
magnetoresistance (AMR) and the planar Hall effect (PHE) are sensitive to the in-plane orienta- 
tion of current and magnetization with respect to crystal axes in a way consistent with the cubic 
symmetry of the system. We also show that the AMR exhibit sign reversal as a function of tem- 
perature, and that it shows significant field dependence without saturation up to 9 T. Our results 
provide a unified description of the anisotropic magnetoresistance effects in epitaxial magnetite 
films and illustrate the need for a full determination of the resistivity tensor in crystalline systems. 
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I. INTRODUCTION 



Magnetite (Fe304) is a 3d transition-metal oxide and an itinerant ferrimagnet with a high 
Curie temperature (858 K) and intriguing transport properties, including a metal insulator 
transition at T ^ 120 K (the Verwey transition) and high spin polarized current attributed 
to its being half metallic. The latter property has given rise to a renewed interest in this 
material as a promising material system for spintronics applications and special efforts 
have been devoted to elucidate its magnet otransport properties. 

Two of the most fundamental magnetotransport properties of itinerant ferromagnets are 
the dependence of the longitudinal (piong) and transverse resistivity (ptrans) on the relative 
orientation of the magnetization (M) and electric current (J), known as anisotropic mag- 
netoresistance (AMR) and planar Hall effect (PHE), respectively. Both phenomena are 
attributed to spin-orbit interaction which mixes spin-up and spin-down states js^. Phe- 
nomenologically, AMR and PHE are commonly expressed as [i^]: 

Plong = P± + (P|| - P±) COS^ ^ (1) 
Ptrans = (P|| " P±) siu (f COS (f, (2) 

where ip is the angle between J and M, and py and p± are the longitudinal and transverse 
resistivities corresponding to M || J and M _L J, respectively. While Eqs. [1] and [2] are 
applicable for isotropic or polycrystalline materials, they fail to fully account for the magne- 



totransport properties of crystalline systems, where additional terms arise [5|]. This has been 
demonstrated for [001]-oriented epitaxial manganites thin films ^, jI, where a description of 
AMR and PHE is provided by taking into account the cubic symmetry of the system. The 
applicable equations are: 

Piong = A cos(2a -2e) + B cos(2a + 29) + C cos(4a) + D (3) 
Ptrans = ^ sin(2a - 29) - 5sin(2a + 29), (4) 

where a and 9 are the angles between M and J relative to the [100] crystal direction, 
respectively (see inset of Fig. [1]). 

In this work we focus our attention on the anisotropic magnetoresistance behaviour of 
epitaxial magnetite films, where the need for an additional term with fourfold symmetry has 
been suggested in the literature ^-lot. We show, however, that Eqs. [3] and [4] are required to 



fully describe the AMR and PHE of cubic epitaxial films for arbitrary directions of current 
and magnetization. 

We find different temperature and magnetic field dependence of the parameters A and 
S, suggesting that they originate from diflFerent scattering mechanisms. The parameter D 
is found to exhibit a strong magnetic field dependence, which is ascribed mainly to the 
eflFect of antiphase boundaries in the Fe304 film. The full description of the AMR and 
PHE behavior of magnetite films allows new insight on the transport mechanisms in this 
compound and enables the optimization of the planar Hall voltage response, which could be 



useful for non- volatile memory devices based on the planar Hall eflFect 
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II. EXPERIMENT 

For this study, Fe304 films ^ 20 nm thick were grown on MgO(OOl) single crystals by 
molecular beam epitaxy in an ultrahigh vacuum deposition system with a base pressure 
of 1 X 10~^ mbar. Prior to film growth the MgO substrate was cleaned under an oxygen 
plasma at 520 K for 30-60 min, which renders the surface well ordered and free of impurities, 
as determined by surface electron diffraction (low energy electron diffraction, LEED, and 
refiection high energy electron diffraction, RHEED) and x-ray spectroscopy (Auger electron 
spectroscopy, AES, and x-xay photoemission spectroscopy, XPS), respectively. For the mag- 
netite growth, an atomic Fe beam was thermally generated from an effusion cell under an 
O2 partial pressure of 2 x 10"'' mbar, with the substrate temperature held at 520 K. Film 
thickness was estimated from a calibrated thickness monitor, while film crystallinity was 
monitored during growth using RHEED. After growth, the film was characterized in situ by 
^^IHEED and XPS, showing that the magnetite films are single crystalline and stoichiometric 
ij]. The MgO substrate crystal orientation was confirmed by Laue diffraction. The films 
were subsequently patterned by photolithography for longitudinal and transverse resistivity 
measurements along different angles (9) of the current with respect to the magnetite crystal 



axes. 
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FIG. 1: Resistivity (p), in log scale, as a function of temperature. Inset: Sketch of the relative 
orientations of the current density J, magnetization M, and the crystallographic axes, piong is 
measured between A and B, ptrans is measured between B and C. 

III. RESULTS AND DISCUSSION 



The temperature dependence of the resistivity is presented in Fig. [T] in semi-log scale. 
We note a change in the slope of the resistivity when approaching the Verwey transition at 
T^ 120 K, indicative of the high quality of the films. Thin magnetite films tend to have 
a less pronounced resistivity discontinuity at the Verwey transition as compared with bulk 
magnetite consistent with earlier transport studies of Fe304 Q. The data in this graph 
essentially corresponds to the temperature dependence of the D term in Eq. [3l 

The variation of piong and ptrans as a function of the angular direction of the applied 
magnetic field (a, relative to [100]Fe3O4) was measured for six different device structures 
with varying angle 9 between the current density direction and the [100] direction of the 
Fe304 film. The measurements were performed with a constant magnetic field of up to 
9 T rotating; in the plane of the film. This field is much larger than the reported in-plane 
magnetization saturating field p^,ll7||; hence, the magnetization is expected to be practically 
parallel to the applied field. Two representative examples of the angular variation of piong 
and ptrans are shown in Fig. [2] for T = 150 K. We find that while ptrans varies sinusoidally, its 
amplitude depends on 9 and the location of the extremal points depend on a, inconsistent 
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FIG. 2: piong and ptrans vs a, the angle between the magnetization M and [100], for different angles 
9 (the angle between the current direction J and [100]) at T=150 K with an apphed magnetic field 
of 9 T for two different samples. The sohd hnes are fits to Eqs. OandHl 

with Eq. O according to which the amplitude is ^-independent. Furthermore, piong lacks 
twofold symmetry. These observations clearly show that Eqs. [H and [2] fail to describe AMR 
and PHE in epitaxial films of magnetite. On the other hand, Eqs. [3] and HJ which take 
into account the crystal symmetry of magnetite, provide a good description for the angular 
dependence of piong and ptrans, as shown by the fits to the experimental data. The equations 
have three parameters and in the fitting process we use the same values for fitting the AMR 
and PHE curves for the six current directions (see Fig. [2j). We note that Eqs. [3] and |4] 
provide a good fit to the data even below the Verwey temperature, where bulk magnetite is 
monoclinic. This is expected since our films are epitaxial and they maintain cubic symmetry 
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FIG. 3: The parameters in Eqs. [3] and [H normalized by the parameter I? as a function of temper- 
ature with an external field of 4 T. The lines are guides to the eye. 

in the plane of the magnetization and current. 

The temperature dependence of the three fitting parameters, 5, and (7, is shown 
in Fig. [31 We note that at high temperatures B is dominant, while at low temperatures 
A and C dominate. The fact that the coefficients A and B show distinct temperature 
dependencies indicate that they have different origins. The terms with the parameter A, 
as their coefficient are identical to the terms that appear in Eqs. [1] and O and they are 
independent of crystal symmetry. On the other hand, the terms with the parameters B and 
C as their coefficients are sensitive to the crystal symmetry. The parameter C appears only 
in the equation for piong and is responsible for a fourfold angular symmetry which provides a 
clear manifestation of crystal symmetry effects. The parameter B contributes to both piong 
and ptrans with a twofold symmetry, and the crystal symmetry is manifested when AMR and 
PHE are measured for current directions away from the principal axes. This may explain 
why the crystal contribution to AMR and PHE at high temperatures has not been reported 
before js]. We note also that, according to Eqs. 3 and 4, the anisotropic magnetoresistance as 
usually defined, AMR = p(M || J) -p(M ± J), yields AMR = 2^ + 25 cos 40, showing that 
for S > ^4 the AMR can be positive or negative, depending of the direction of the electric 
current (9). In particular, for = 7r/4, AMR = 2{A — B) ^ and the AMR, as conventionally 
defined, changes sign at around 150 K. 
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The field dependence of the magnetoresistance fitting parameters was also measured at 
selected temperatures and representative results for one sample are shown in Fig. [31(a). 
While the field dependence of A and B is weak, the parameter C varies strongly with 
the magnetic field amplitude, with no sign of saturating up to 9 T. The field dependence 
yields quite dramatic changes in the AMR as can be seen in Fig. [11(b), which exhibits a 
much stronger fourfold symmetry component when the field is increased from 2 T to 9 T. 
The field dependence of the magnetoresistance parameters varies with temperature, but the 
qualitative behavior is similar to that shown in Fig. [11(a). 

The parameter which represents magnet oresistivity term insensitive to a and 0, is also 
found to vary strongly with the magnetic field, giving rise to large magnetoresistances of 
the order of -10% at 9 T. While the magnetoresistance ratio of bulk magnetite is relatively 
weak, of the order of ^0.1% fl5^, much larger eflFects, of up to -10% at 5 T jist, are observed 
in thin films. 

Several different factors may aflFect the magnetotransport behavior of magnetite, including 
the spin-orbit coupling, electron localization on approac hing the Verwey transition, changes 
in the spin polarization as a function of temperature |19|, and changes in the character 
of the charge carrier at temperatures near room temperature [20|. The combination of all 
these factors makes a correct identification of the dominant scattering processes in play 
challenging. We do note, however, that the difference between bulk compounds and thin 



films grown on MgO (particulary, lack of 
presence of antiphase boundaries (APB) 



ligh field saturation) has been attributed to the 
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23| . The APB disturb the magnetic 



configuration even at high magnetic fields, particularly close to the boundary, as illustrated 
in the inset of Fig. [51(b); hence, they are likely to be the origin of the magnetic field 
dependence of the parameter D. We turn now to examine whether there could be a link 
between APB and the observed field dependance of AMR. 

We follow previous treatments and calculate the spin configuration, /3(x), in an infinite 
ferromagnetically-coupled spin chain with anti-ferromagnetic boundary at the origin. The 



energy per unit area of such a spin chain to the left to the origin is given by [24 1 

s= r [-MsHco8f3 + AF(^f]dx (5) 
J-oc dx 

where the first term in the integral is the Zeeman energy and the second term is the near- 
est neighbor exchange contribution. The angle between the magnetization, M^, and the 
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FIG. 4: (a): (right axis) The parameter D and (left axis) the other parameters from Eqs. [3]and[ll 
normahzed by D as a function of the external field at T=150 K. The lines are guides to the eye. 
(b): AMR normalized to its maximum value as a function of a for two different external magnetic 
fields at T = 150 K. The solid hues are fits to Eq. [3 

magnetic field, i7, is denoted by /3, and is the exchange stiffness constant. We neglect 
the magnet ocrystalline, since it is much smaller than the external field that we used in our 
measurements p^, 117|. Standard variational calculus procedures yield: 

/3(x) =4arctan(Ce^/^°) (6) 

where = 2Af/MsH^ C = tan(/3^i;^/4) and 2/3af is the angle between the spins on both 




FIG. 5: The field dependence of a2 (the numerical average of cos 2/3) and a4 (the numerical average 
of cos 4/3) for the whole domain and (Inset) for the vicinity of the boundary alone (see text), (b): 
an illustration of spin orientation of a ferromagnetic chain with anti-ferromagnetic coupling at the 
origin, taken from [l8]. 

sides of the boundary which depends on the anti-ferromagnetic coupling, and can be calcu- 
lated based on the equilibrium condition of the chain |l8|. We note that this model does 
not take into account the distribution in the orientation of the APB, since /3(x) does not 
depend on the angle between the APB and magnetic field. 

Figure [5] presents the spatial average of cos 2/3 and cos 4/3 denoted as (cos 2/3) and (cos 4/3), 
res pec tively, using a ferromagnetic coupling of 3 K and anti-ferromagnetic coupling of -22 



K |25|. The averages are as = 5^ cos(2/3(jrf)) and = ^;^ cos(4/3(jrf)) 

respectively, where d is the separation between spins (about 1 A j22|) and the domain size 
is ^ 350 A (26! . Changing the domain size in a reasonable range between 250 A and 350 A 
does not affect the results significantly. 



The field dependence of a2 anc 
C. However, as suggested before 



Oa is much smaller than the field dependance of B and 



181], the largest contribution to resistivity is close to the 
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boundary (due to the significant change in the spin orientation) and taking the average only 
in the vicinity of the boundary yields a more dramatic field dependance (see inset of Fig. E]). 
Thus, one can not exclude that a more sophisticated model that would consider the effect 
of spin variation on p would yield better agreement with experiment. 

In summary, we measured the temperature and magnetic field dependence of the longi- 
tudinal and transverse resistivities of Fe304 as a function of the current direction relative 
to the crystal axes. The results shed light on the interplay between magnetism and elec- 
trical transport in this class of materials and may serve as a basis for further study of the 
microscopic origin of magnetotransport properties of magnetite. 
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